function ordertest()

  X_GRID_SIZE = 32;
  Z_GRID_SIZE = 32;
  L = 100;
  TIME_FINAL = 100;
  TIME_STEPS = 200;
  fm = [1:2]; % frequency multiplier

  pwrs = [0:0]
  pwrs = 2.^pwrs;
  amp = 0.01 .* pwrs

  evalpoints = [10:5:95];


  metricorder = [0 1 2];


  %res0 = zeros(length(fm), length(amp), length(evalpoints));
  res1 = zeros(length(fm), length(amp), length(metricorder), length(evalpoints), );

  for fmi=1:length(fm)
    for ampi=1:length(amp)
      disp(sprintf('*** freq: %d,  amp: %d', fm(fmi), amp(ampi)));
      r = h2(X_GRID_SIZE, L, Z_GRID_SIZE, TIME_FINAL, TIME_STEPS, amp(g), fm(z), evalpoints, metricorder);
      res0(z, g, :) = r;
      disp(r)
    end
  end

  save('ordertest.mat', 'res0', 'resy', 'fm', 'amp');





end
